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基于 叶绿体 基因 证 据 的 民族 药 演 白 珠 复合 群 系统 发 育 关 系 
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摘要 : 汗 白 珠 是 我 国 重 要 的 民族 药 用 植物 ,广泛 分 布 于 长 江 以 南 地 区 ， 是 一 个 分 类 困难 的 复 
合群 ,为 了 探讨 其 种 下 分 类 关系 ,我 们 对 演 白 珠 复合 群 ARBAA, ARARA AER 
等 三 个 变种 ) 进行 网 罗 式 采样 ， 基 本 覆盖 了 该 复合 群 在 中 国 的 分 布 范 围 ,同时 包括 菲律宾 和 
马来西亚 的 各 1 个 居 群 ， 共 计 81 个 居 群 241 个 个 体 ， 通 过 联合 两 个 变异 位 点 适中 的 叶绿体 
片段 rp133-psaJ 和 trnL-rp132， 构 建 基 于 最 大 似 然 法 (Maximum Likelihood) 和 贝 叶 斯 法 
(Bayesian Inference) 的 系统 发 育 树 ， 以 及 Neighbor-Net 法 构建 系统 发 育 网 络 。 研 究 表 明 : 
演 白 珠 复合 群 内 具有 明显 的 遗传 差异 ,这 种 差异 性 相 比 形态 , 与 地 域 分 布 相关 性 更 大 。 系 统 
发 育 分 析 显 示 : 溃 白 珠 复合 群 分 为 三 个 支 系 , 一 文 包括 台湾 和 菲律宾 的 居 群 ， 为 变种 秀 果 
珠 ， 符 合 前 人 分 类 结果 ; 一 文 包 括 分 布 于 横断 山脉 区 域 的 居 群 ， 由 变种 毛 演 白 珠 和 演 白 珠 组 
成 ; 剩余 一 文 包 括 分 布 于 华东 南 区 域 的 居 群 ， 同 样 由 变种 毛 汗 白 珠 和 污 和 白 珠 组 成 。 分 子 证 据 
文 持 基于 形态 分 类 的 秃 果 白 珠 作为 变种 ， 而 不 支持 毛 滇 白 珠 作为 变种 的 处 理 。 这 样 的 遗传 分 
化 式样 可 能 是 由 地 理 隔离 导致 , 这 一 结果 为 民族 药 滇 白 珠 的 种 质 资源 评价 和 保护 利用 奠定 了 
分 类 学 框架 。 然而 ， 本 研究 变种 溃 白 珠 居 群 间 的 谱系 发 生 关系 尚未 得 到 很 好 地 解决 ， 特别 是 
无 法 明晰 华东 南大 部 分 居 群 间 的 遗传 分 化 格局 ， 需 要 进一步 补充 样本 通过 谱系 地 理学 手段 ， 
并 寻找 分 辩 率 更 高 的 核 基因 分 子 标记 进行 更 深入 的 研究 。 
关键 词 : 溃 白 珠 复合 群 ， 叶 绿 体 DNA， 系 统 发 育 ， 种 内 遗传 差异 ， 种 内 界定 ， 地 理 隔离 
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Abstract: As an important ethnical medicinal plant in China, Gaultheria leucocarpa var. 
yunnanensis is a taxonomically difficult complex and widely distributed throughout south of the 
Yangtze River. In order to understand its intraspecific relationship, we sampled 241 individuals in 
81 populations of G. leucocarpa var. yunnanensis complex (1.e., G. leucocarpa var. crenulata, G. 
leucocarpa var. psilocarpa and G. leucocarpa var. yunnanensis) representing all Chinese varieties 
and covering almost recorded distribution ranges in China (Mainland and Taiwan). One population 
of G. leucocarpa var. psilocarpa from the Philippines and one population of G. /leucocarpa var. 
leucocarpa from Malaysia were also included. Based on the data from two plastid DNA loci with 
variable sites in appropriate number, i.e., rpi33-psaJ and trnL-rpl32, we reconstructed 
phylogenetic trees using both Maximum Likelihood and Bayesian Inference methods, and 
phylogenetic network using Neighbor-Net method. Our results showed that there was a distinct 
genetic differentiation in the G. leucocarpa var. yunnanensis complex, and this differentiation was 
more correlated with geographical distribution rather than morphology. This complex is mainly 
divided into three clades: one clade is composed of populations from Taiwan and the Philippines 
that are previously categorized as G. /leucocarpa var. psilocarpa; one clade includes populations 
of the G. leucocarpa var. yunnanensis and G. leucocarpa var. crenulata with a distribution in the 
Hengduan Mountains region; and the remaining clade includes the populations mainly distributed 
in the southeastern range of China, which is composed of G. leucocarpa var. yunnanensis and G. 
leucocarpa var. crenulata as well. Molecular evidence therefore supports G. leucocarpa var. 
psilocarpa was treated as a variety but does not support G. /leucocarpa var. crenulata as a variety. 
Such genetic differentiation pattern may be caused by geographical isolation. This result lays a 
basic taxonomic framework for the evaluation and protection of germplasm resources for this 
ethnical medicinal plant. However, the relationships among most mainland populations of G. 
leucocarpa var. yunnanensis has not been well resolved in this study. In particular, it is still unable 
to clarify genetic differentiation pattern among most populations in southeastern China. Therefore, 
it is necessary to further carry out phylogeographic study at population level using higher variable 


nuclear markers in future. 


Key words: Gaultheria leucocarpa Bl. var. yunnanensis complex, cpDNA, phylogeny, 
intraspecific genetic difference, intraspecific delimitation, geographic isolation 
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白 珠 Gaultheria leucocarpa Bl. var. yunnanensis (Franch.) T. Z. Hsu & R. C. Fang 是 白果 
白 珠 G. leucocarpa W]2E Th, RJE T 4EBS 46 FL Ericaceae 越 桔 亚 科 Vaccinioideae 白 珠 树 族 
Gaultherieae 白 珠 树 属 Gaultheria Kalm & L. (Kron et al., 2002) 。 在 国内 广泛 分 布 于 长 江 流 
域 以 南 各 省 区 、 包 括 台 湾 和 海南 岛 ， 生 长 在 海拔 3500 米 以 下 的 针 阔 混交 林 下 。 演 白 珠 是 我 
国 重要 的 民族 药 用 植物 ， 为 很 多 地 区 特别 是 西南 地 区 ) WARR RR IR RER 
族 、 纳 西 族 、 水 族 、 瑶 族 、 壮 族 、 布 依 族 等 10 个 民族 习 用 〈 马 小 军 等 ，2001) 。 也 是 我 国 
医药 界 的 战略 性 资源 ， 市 场 需求 广泛 ， 该 植物 含 多 种 化 学 成 分 ， 目 前 被 开发 的 相关 药物 有 主 
眩 尝 症 的 ' 溃 白 珠 糖浆 "(批准 文 号 :; 国药 准 字 Z20040088) ， 治 疗 风 湿 肿 痛 的 ' 金 骨 血 片 ， 
(批准 文 号 : 国药 准 字 Z20080203) ， 具 有 活血 祛 病 、 消 肿 止 痛 功 效 的 复方 透 骨 香 软膏 '( 批 
准 文 号 : 国药 准 字 Z20026850) 等 ; 并且 其 枝 、 叶 富 含 芳香 油 (水 杨 酸 甲 酯 ， 达 0.5-0.8%) ， 
是 工业 化 提取 精油 的 良好 原料 。 

与 汗 白 珠 同属 于 白 珠 树 属 的 芳香 白 珠 G. fragrantissima Wall., 五 雄 白 珠 G. semi-infera (C. 
B. Clarke) Airy Shaw. £L HK G. hookeri C. B. Clarke、 尾 叶 和 白 珠 G. griffithiana Wight 等 ， 
因 其 不 同 程度 地 含有 水 样 酸 甲 酯 而 被 混 作 演 白 珠 的 基 原 使 用 〈 乔 文 林 ，2013) ， 但 是 这 些 植 
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物 均 不 售 重 要 抗 炎 活性 成 分 一 溃 白 珠 树 苷 (Gaultherin) (Liu et al., 2013) 。 上 述 四 种 植物 
与 汗 白 珠 在 形态 上 具有 明显 差异 ， 容 易 辨 识 ， 药 源 鉴 定 不 会 造成 较 大 困难 。 然 而 , 溃 白 珠 基 
原本 身 的 次 生 代 谢 产物 种 类 及 含量 差异 性 较 大 ， 赵 玉 娟 等 (2002) 、 马 小 军 等 (2002) 、 乔 
文 林 (2013) ~ Liu etal. (2015) 等 对 不 同 基 原 / 产 地 的 滇 白 珠 开展 了 化 学 成 分 相关 研究 ， 发 
现 化 合 物 及 其 含量 具有 显著 地 域 性 差异 ， 为 该 植物 的 资源 开发 利用 提供 了 一 定 的 理论 参考 。 
但 是 ,这 些 工 作 主 要 聚焦 在 药物 化 学 成 分 含量 的 研究 上 ,往往 忽略 了 基 原 调查 ， 以 及 生物 学 
(包括 遗传 学 ) 特性 与 地 理 分 布 关系 的 研究 。 
就 溃 白 珠 本 身 而 言 ， 该 植物 形态 变异 不 大 但 遗传 分 化 显著 〈 陆 露 ，2009) ， 是 一 个 分 类 
困难 的 复合 群 ,是 水 平和 垂直 空间 分 布 最 广 的 国产 白 珠 树 属 植 物 , 其 余 的 种 均 集 中 分 布 于 喜 
马 拉 雅 -横断 山区 (Fang & Stevens, 2005) 。 除 汗 白 珠 外 ， 和 白果 和 白 珠 还 包括 其 它 5 个 变种 ， 
我 国境 内 有 毛 演 白 珠 G. leucocarpa var. crenulata (Kurz) T. Z. Hsu 和 秃 果 白 珠 G. leucocarpa 
var. psilocarpa (Copeland) R. C. Fang (DATE) 2 个 变种 (Fang & Stevens, 2005) 。 这 
3 个 国产 变种 果实 均 为 黑 紫 色 ， 构 成 了 污 和 白 珠 复合 群 ， 均 可 入 药 ， 且 均 为 民族 药 滇 白 珠 的 基 
原 。 汗 和 白 珠 复合 群 的 叶 型 、 毛 被 、 花 部 结构 和 果实 形态 有 一 定 程 度 变异 但 不 显著 ， 其 分 类 处 
理 一 直 以 来 倍 受 争议 (Copeland，1932; Sleumer, 1967; Ying, 1976; 徐 廷 志 ，1991; Fang 
& Stevens, 2005; Fritsch etal., 2008) 。 除 毛 演 白 珠 明显 被 毛 外 ， 其 余 2 个 国产 变种 ， 即 
溃 白 珠 、 秀 果 白 珠 在 形态 上 较 难 相互 辨别 ， 仅 依靠 萌 果 是 否 被 毛 所 区 别 。 同 时 ， 汗 白 珠 变种 
内 形态 变异 有 时 大 于 变种 间 变 异 ， 为 其 分 类 学 处 理 增加 了 难度 。 

演 白 珠 复合 群 有 具有 明显 的 遗传 分 化 ， 广 域 分 布 的 滇 白 球 和 中 国 植物 志 英 文 版 (Fang & 
Stevens, 2005) 中 作为 汗 白 珠 异 名 处 理 的 屏 边 白 珠 G. leucocarpa var. pingbienensis C. Y. Wu 
ex T. Z. Hsu〈 狭 域 分 布 于 云南 省 屏 边 县 ) 形态 相近 。 近 年 来 ,利用 分 子 手段 对 滇 白 珠 复合 群 
展开 的 研究 揭示 了 复合 群 内 部 分 变种 的 系统 关系 : 叶绿体 基因 trnL-trnF 和 rpll6 序列 中 的 一 
些 插 入 /缺失 对 谱系 界定 起 到 了 重要 作用 ， 联 合 分 析 高 度 支 持 了 变种 屏 边 白 珠 与 白果 白 珠 原 
变种 构成 一 支 , 而 变种 溃 白 珠 却 与 形态 有 别 的 变种 毛 汗 白 珠 聚 成 男 一 支 , 这 两 个 分 支 构成 姐 
妹 群 关系 〈 陆 露 ，2009) 。 这 些 结果 显示 出 滇 白 珠 复合 群 形态 和 遗传 的 不 一 致 ， 基 于 传统 的 
形态 分 类 无 法 明晰 滇 白 珠 复 合群 的 种 内 分 类 问题 。 经 野外 调查 发 现 , 汗 白 珠 和 毛 滇 白 珠 虽 然 
在 形态 上 以 毛 被 区 分 , 但 两 者 分 布 范 围 和 生长 环境 相似 ,常常 混 生 ,而 且 一 些 溃 白 珠 的 幼 株 
茎 时 也 常常 被 毛 。 因 此 ， 了 解 滇 白 珠 复合 群 的 遗传 分 化 格局 对 重新 认识 基于 形态 划分 的 变种 
有 具有 重要 意义 ， 叹 需 对 该 类 群 开展 遗传 分 化 格局 与 地 域 、 形 态 的 相关 性 研究 。 本 研究 首次 对 
溃 白 珠 复合 群 在 中 国境 内 进行 了 大 范围 取样 ， 几 乎 涵盖 了 中 国境 内 演 白 珠 复合 群 的 分 布 区 
域 ， 同 时 包括 东南 半岛 的 菲律宾 和 马来西亚 各 1 个 居 群 ， 共 计 81 个 居 群 Qd, Bp. B 
中 ， 在 云南 武 定 (WD) 和 贵州 印 江 (YHA) 分 别 采集 各 1 个 毛 溃 白 珠 居 群 ， 台 湾 (HXQ) 
MIER F) 分 别 采 集 各 1 个 秀 果 白 珠 居 群 ， 以 及 马来西亚 (WH) 的 1 个 白果 白 珠 原 变 
种 居 群 。 选用 进化 速率 相对 快 的 叶绿体 DNA. 非 编码 区 序列 rp133-psaJ 和 trnL-rpI32. 首次 开 
展 汗 白 珠 复 合群 的 遗传 分 化 和 种 下 系统 发 育 关系 研究 , 进而 为 汗 白 珠 基 原 界定 和 之 后 化 学 活 
性 成 分 的 地 域 性 研究 提供 遗传 学 理论 框架 。 


= 
p= = 
— 

9 
C) 


1 材料 和 方法 
1.1 实验 材料 

本 研究 共 收 集 了 汀 白 珠 复合 群 81 个 居 群 和 1 个 芳香 白 珠 居 群 〈 作 为 外 类 群 ) 的 分 子 材 
料 ， 分 别 采 自 云南 (17， 其 中 WD 居 群 为 毛 汗 白 珠 ) 、 贵 州 (20， 其 中 YHA 居 群 为 毛 汗 白 
珠 ) 、 四 川 (G5 、 重 庆 COD 、 湖 北 G) WA (11) 、 江 西 (4) 、 福 建 (4) 、 广 东 (9)、 
广西 (5) 、 台 湾 〈1， 秀 果 白 珠 ) 、 菲 律 宾 〈1， 秃 果 白 珠 ) 、 马 来 西亚 〈1， 白 果 白 珠 原 变 
种 )， 有 76 个 汗 和 白 珠 居 群 ; 芳香 白 珠 (FX) 采 自 云南 大 理 。 每 个 居 群 随机 挑选 3 个 个 体 ( 注 : 
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峨眉 山 EMS) 仅 发 现 一 个 个 体 ) ， 包 括 外 类 群 共计 244 个 个 体 〈 表 1) 。 采 和 集 材 料 时 ， 居 
群 内 每 个 个 体 至 少 相 隔 10 米 以 上 。 新 鲜 叶 片 立 即 用 变色 硅胶 干燥 ， 和 凭证 标本 储存 于 昆明 医 
科大 学 药学 院 标本 储藏 室 。 
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Table 1 Sampling and voucher information of populations of the Gaultheria leucocarpa var. 


变种 名 称 


Name of varieties 


AAR 


Gaultheria 


yunnanensis complex and related taxa 


leucocarpa 


var. yunnanensis 


居 群 编号 采集 地 采集 号 个 体 数 
Population L : Collection Samples 
ocation 
code number (n) 
国 云南 安宁 
AN d . d . LRF-LU1001 3 
Anning, Yunnan, China 
中 国 云南 保山 昌 字 
BS | B LJ-Liul001 3 
Changning, Baoshan, Yunnan, China 
国 四 川 布 托 
BT T [783 . LL-2017-02 3 
Butuo, Sichuan, China 
中 国 贵 州长 顺 
CS li | | LL-2017-19 3 
Changshun, Guizhou, China 
中 国 云南 大 理 苍 ， 
DL . LL-2019-03 3 
Cangshan, Dali, Yunnan, China 
屏 边 大 围 山 
DWS iioii THU f LRF, LU002 3 
Daweishan, Pingbian, Yunnan, China 
国 贵 州都 匀 
DY13 TEST i LL-2017-13 3 
Duyun, Guizhou, China 
贵州 都 们 
DY14 TEE aA l LL-2017-14 3 
Duyun, Guizhou, China 
国 广西 大 瑶 山 
DYS ü ^ : . LL-2012-03 3 
Dayaoshan, Guangxi, China 
到 四 川 峨眉 山 
EMS RITU | iin f LL-2017-05 1 
Emeishan, Sichuan, China 
国 云南 富源 
FY a LL-2016-23 3 
Fuyuan, Yunnan, China 
中 国 贵州 兴 仁 
GDH a l LL-2017-21 3 
Xingren, Guizhou, China 
中 国 贵州 贵 定 
GD ed | LL-2017-16 3 
Guiding, Guizhou, China 
中 国 贵州 龙 里 
GLL EM | LL-2017-17 3 
Longli, Guizhou, China 
中 国 贵州 贵阳 
GY T LL-2017-18 3 
Guiyang, Guizhou, China 
中 国 贵州 珠 市 
GZS | LL-2016-20 3 
Zhushi, Guizhou, China 
中 国 贵 州 赫 章 
HZ i i | LL-2016-14 3 
Hezhang, Guizhou, China 
中 国 云南 景 东 
JD : LL-0609 3 
Jingdong, Yunnan, China 
JWS 中 国 广 西柳 州 九 万 山 LL-2012-01 3 
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1.2 引物 的 筛选 


在 引物 筛选 阶段 , 我 们 从 白果 白 珠 变种 已 有 各 居 群 中 根据 地 理 相对 距离 较 远 原则 挑选 出 


ZABI CRA). SARE BAAR HIPPIE GRAO AREL CRARO, 
台湾 玉山 ( 秀 果 白 珠 ) 、 马 来 西亚 云顶 (白果 白 珠 原 变种 〉 共 6 个 居 群 ， 每 个 居 群 随机 选取 


2 个 个 体 ， 共 12 个 个 体 。 从 杜 鹏 花 科 近 缘 的 山茶 科 Theaceae 中 ,选取 了 进化 速率 较 快 的 10 


个 叶绿体 基因 片段 (P37、P45、P57、P58、P60、P61、 


P64, P69, P70 (HU: 7p133-psaJ) 、 


P72 , Yang et al, 2013 ) 进行 分 析 ， 最 终 选 出 7p133-psaJ C F: 5 
TCCGGATGCGTTAACATTCCCCTT 3'; R: 5' CCTTGGAAGGGTAACACACAGGTGC 3^) 和 


在 杜 鹏 花 科 中 拥有 较 高 分 辨 率 的 trnL-rp132 (trnL:5 
rpL32:5' CAGTTCCAAAAAAACGTACTTC 3', Shaw et 
分 子 片段 。 


1.3 DNA 提取 与 PCR 扩 增 
植物 总 DNA 的 提取 使 用 改进 的 4xCTAB 法 ， 之 后 


CTGCTTCCTAAGAGCAGCGT 3* 
al.，2007) 作为 开展 系统 发 育 研 究 的 


J 1% 的 琼脂 糖 凝 胶 电 泳 检测 ， 根 据 


DNA 主 带 的 整齐 度 和 亮度 判断 DNA 的 质量 , 继而 用 紫外 分 光 光 度 计 (NanoDrop@ ND-1000, 
Thermo Fisher Scientific, Delaware, USA ) 测 定 DNA 的 浓度 和 纯度 , 后 用 ddH2O 稀释 至 30-50 
nge uL! 以 备 后 续 PCR 扩 增 使 用 。 扩 增 体 系 为 : 10 pL PCR Master Mix (2x) (Thermo 
Scientific) ，9 uL ddH2O， 正 反 向 引物 C10 ng*uL % 0.4 uL. 1 kL 模板 DNA。PCR 扩 
增 反 应 程序 为 : 94 C 预 变性 4min; 94 C 变性 45s, 52 C CrnL-rpl32)/60 C CrpI33-psaJ) 


退火 1 min, 72 'C 延伸 45 s, 循环 35 次; 72 C. 延 人 
琼脂 糖 凝 胶 电泳 ,检测 扩 增产 物 的 产量 和 质量 。 剩 余 的 


10 min。 取 2 uL 扩 增 产物 进行 1.5% 
PCR 扩 增 产物 纯化 后 , 采用 5.5 pL 


的 反应 体积 ， 包 括 1.5 pL 的 测序 混合 物 (BigDye )，1 kL 的 测序 引物 ，1 uL 的 纯化 PCR 产 
物 和 2 uL 的 双 节 水 。 反 应 程序 为 : 95 'C 预 变性 40s; 95 'C 20s. 50 'C 10s、60 'C 1 min 30 
s, d 35 个 循环 进行 测序 反应 。 反 应 产物 经 过 75% 乙 醇 、 无 水 乙醇 纯化 后 ， 加 入 10 uL 的 
双 蒸 水 充分 溶解 ， 用 ABI 3730 xl 自动 测序 仪 CApplied Biosystems, Foster City, California, 


USA) 进行 测序 。 


1.4 数据 分 析 


所 得 序列 用 Sequencher 软件 进行 拼接 。 拼 接 好 的 序列 首先 逐一 对 每 一 个 变异 位 点 的 原 
始 峰 图 进行 校对 ， 确 认 无 误 后 用 Se-al 软件 进行 比 对 构建 矩阵 ， 根 据 头 尾 的 缺失 情况 对 序列 
做 适当 调整 。 拢 阵 中 的 插入 /缺失 〈indel) 处 理 为 单 次 进化 事件 。 为 解析 目前 所 获 各 居 群 间 


的 系统 进化 关系 ， 我 们 对 两 个 叶绿体 片段 rpI33-psaJ 和 


trnL-rpl32 进行 了 联合 分 析 。 本 研究 


用 最 大 似 然 法 (Maximum Likelihood, ML) 和 贝 叶 斯 法 (Bayesian inference, BI) 构建 系统 


树 , 以 比较 不 同 拓扑 结构 的 异同 , 并 计算 支持 率 (Stamatakis et al, 2008). 利用 软件 jmodeltest 
2.110 (Darriba, 20120, 分 别 检 测 到 两 个 叶绿体 基因 片段 的 最 佳 碱 基 蔡 代 模 型 都 为 GTR。 最 
大 似 然 法 依据 DNA 序列 突变 是 随机 的 这 种 假设 ,将 构建 系统 树 所 用 到 的 分 枝 长 度 与 推断 肯 
定 发 生 的 额外 突变 结合 起 来 ， 利 用 特定 突变 概率 的 所 有 假设 共同 组 成 一 个 演化 模型 (model 
of evolution )， 选 取 具 有 最 大 可 能 性 的 树 作为 引证 ， 分 析 其 系统 进化 关系 ， 并 计算 支持 率 


(Bootstrap ，BP)。 贝 叶 斯 法 构建 系统 进化 关系 树 时 ， 


A i 


"Hog E 


需 对 4 个 平行 进化 链 进 行 1000 万 代 


算 , 每 10000 代 保 存 一 棵 进化 树 , 分 别 去 掉 最 初 和 最 后 的 25 % 系 统 进化 树 作为 “burn-in”。 
个 分 支 的 支持 率 由 保存 的 进化 树 中 出 现 的 后 验 概率 (Posterior probability; PP) 表示 。 后 
概率 值 PP < 0.95 视 为 低 或 较 弱 的 支持 ，0.95 一 1.00 视 为 高 或 好 的 支持 。 两 种 建树 方法 都 
过 CIPRES 网 上 建树 平台 (https://www.phylo.org/) 实现 ， 最 后 得 到 的 基于 ML 和 BI 法 的 
格 一 致 树 用 Figtree 1.43 软件 Chttp;//tree.bio.ed.ac.uk/software/figtree/) 查看。 为 了 直观 地 


判断 滇 白 珠 复 合群 的 系统 演化 关系 ， 用 Splitstreed 软件 运用 邻接 网 法 (Neighbor-Net) 构建 
系统 发 育 网 络 (Bryant, 2003), Bootstrapping 通过 1 000 次 重复 ， 支 持 率 大 于 50% 的 显示 在 
Es 


2 结果 与 分 析 


本 研究 成 功 地 从 溃 白 珠 复合 群 241 个 样本 和 外 类 群芳 香 白 珠 3 个 样本 获得 了 rp133-psaJ 

和 trnL-rpl32. 片段 .rp133-psaJ 片段 比 对 和 矩阵 长 678 bp, 共有 241 条 序列 , 一 致 位 点 (conserved 

site) 有 665 个， 变异 位 点 (variable site) 有 13 个 ， 含 简约 信息 位 点 (parsimony informative 
site) 有 12 个 ; trnL-rp132 片段 比 对 矩阵 长 599 bp， 共 有 241 条 序列 ， 一 致 位 点 有 592 个 ， 
变异 位 点 有 7 个 ， 含 简约 信息 位 点 有 $ 个 。 基 于 两 个 叶绿体 片段 rp133-psaJ 和 trnL-rpI32 的 
联合 分 析 和 矩阵 长 度 为 1277 bp， 一 致 位 点 有 1257 个 ,变异 位 点 有 20 个 , 含 简约 信息 位 点 有 

17 个 。 汗 和 白 珠 复合 群 联合 矩阵 的 G+C 含量 为 30.70%~30.90%， 其 中 包括 一 段 在 trnL-rp132 
序列 上 的 、 长 度 为 31 bp 的 插入 /缺失 序列 AATTGAATTTCATATTGATTAGTCCAAAACA; 
该 插入 /缺失 序列 为 台湾 居 群 (HXQ) 所 特有 。 另 外 ，3 条 长 1215 bp 的 芳香 白 珠 序列 ，G+C 
含量 为 31.11% ， 该 序列 亦 存 在 一 段 31 bp 的 插入 /缺失 序列 

AATTGAATTTCATATTGATTAGTCCAAAACA (5 HXQ 居 群 的 插入 /缺失 序列 一 致 但 位 置 
不 同 )。 图 1 显示 了 贝 叶 斯 法 (BI) 构建 的 系统 发 育 树 ， 该 树 与 通过 最 大 似 然 法 (ML) 构建 
的 树 拓扑 结构 一 致 ， 在 BI 树 上 标注 了 ML 树 的 支持 率 ，Neighbor-Net 构建 的 滇 白 珠 复合 群 

- 系统 发 育 网 络 关系 见 图 2. 

Y 基于 ML 和 BI 分 析 可 得 到 以 下 结果 : CIO 广 域 分 布 的 演 白 珠 复合 群 存在 丰富 的 遗传 变 

异 。ML 和 BI 系统 树 的 拓扑 结构 基本 一 至， 汗 白 珠 复合 群 大 致 可 以 分 为 三 个 大 的 支 系 (BP 

二 100%，PP 二 1.00， 图 1)。 台 湾 (HXQO MIRE C) 的 秀 果 白 珠 居 群 在 较 基部 的 位 置 

^ AUS — 4x (A) (BP 王 91%，PP 王 1.00)， 是 较 早 分 化 出 来 的 类 群 ， 与 其 它 所 有 居 群 相 区 
ft 别 ， 随 后 两 个 居 群 内 的 个 体 又 各 自 聚 为 一 支 (HXQ: BP=100%, PP—1.00; F: BP 二 98%， 
,一 PP=1.00); 剩余 居 群 可 以 分 为 两 个 主要 分 支 ， 第 二 分 支 为 分 布 在 横断 山脉 的 云南 和 四 川 、 

^ 华东 南 的 广西 、 广 东 、 福 建 的 部 分 个 体 以 及 马来西亚 CWE 居 群 聚 为 的 一 支 ， 具 有 较 高 支 

| 持 率 (BP 一 90%，PP 一 1.00)。 这 一 分 支 进而 分 为 4 个 小 的 亚 分 支 ; 马来西亚 的 白果 白 珠 原 
= 变种 (WH) 位 于 较 基部 的 位 置 , 单独 形成 一 支 (B1) 与 其 他 居 群 分 开 (BP 二 96%, PP— 1.00), 

© 广西 大 瑶 山 (DYS-1) HJ RÆ CDYDS-1, DYDS-2) XJj— 3 (B2) (BP= 9796, PP 
— 1.000, 福建 龙岩 (FDQS-1) 和 福建 漳州 (FMHS-1) 聚 为 一 支 (B3) (BP —999$, PP— 1.00), 

剩 下 的 个 体 分 别 来 自 云 南 和 四 川 的 AN (2) DL (1) 、QJ (1) WD G), EMS (D 、 

XC (3) 、SB G) 居 群 ， 为 该 大 文 的 主要 成 分 (B4) 。 其 它 剩 余 的 居 群 聚 为 第 三 分 支 (BP 

—8996, PP—1.000, ix —4) cT 3 个 小 的 亚 分 支 ， 第 一 文 (C1) 是 位 于 较 基部 的 屏 边 

大 围 山 的 两 个 个 体 (DWS-1，DWS-2)， 得 到 了 较 高 的 支持 率 (BP 二 98%，PP 二 1.00), 福建 

龙岩 两 个 居 群 的 两 个 个 体 (FGT-3 和 FHJS-3 ) 聚 为 一 支 (C2) (BP=95%, PP—1.000, RF 

的 居 群 个 体 构 成 一 个 大 的 分 支 (C3) ， 序 列 占 总 数 的 80% 以 上 “201 条 序列 )， 这 些 居 群 个 体 

之 间 的 关系 尚未 得 到 很 好 地 解决 ， 只 有 贵州 石 哇 〈《SWT-1) 和 云南 曲靖 的 〈ZY-2) RIX 

(CA) 而 与 其 他 序列 分 开 (BP—959$, PP—0.97). Neighbor-Net 的 网 状 结构 显示 汗 白 珠 复 合 

群 也 可 以 大 致 分 为 三 个 支 系 (图 2)。 人 台湾、 菲律宾 和 马来西亚 的 居 群 以 及 分 布 在 横断 山脉 

区 域 的 居 群 聚 为 一 支 (红色 圈 注 )， 支持 率 较 低 ; 与 ML 和 BI 系统 发 育 分 析 结 果 一 致 ， 台 湾 
的 秃 果 白 珠 居 群 和 地 理 上 相距 更 近 的 福建 汗 白 珠 居 群 未 聚 在 一 支 。 大 多 数 来 自 华 东南 地 区 的 

居 群 聚 为 一 支 〈 蓝 色 圈 注 )。 汗 东南 地 区 来 自 屏 边 大 围 山 的 居 群 (DWS) 单独 构成 一 支 〈 紫 

BEE). 

(2) 漂白 珠 复合 群 的 遗传 式样 和 毛 被 这 一 表 型 式样 不 相符 。ML 和 BI 分析 均 显 示 采 


202001.00089v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


云南 武 定 (WD) 和 贵州 印 江 CYHAO 的 毛 滇 白 珠 居 和 群 并 没有 聚 为 一 支 而 是 分 别 和 采 自 相同 
区 域 的 汗 白 珠 居 群 构 成 一 支 ( 图 1), Neighbor-Net 结果 也 显示 毛 溃 白 珠 与 大 部 分 相同 区 域 的 
演 白 珠 居 群 聚 在 一 起 《图 2)。 
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图 1 基于 贝 叶 斯 法 (BI) 构建 的 滇 白 珠 复合 群 系统 发 育 树 ， 外 类 群 为 芳香 白 珠 。 系 统 发 育 树 上 标注 了 基于 
ML 分 析 的 支持 率 与 BI 分 析 的 后 验 率 : 大 于 或 等 于 50% 的 支持 率 BP， 分 校 上 左 侧 ) 和 大 于 或 等 于 50% 
的 后 验 率 PP， 分 枝 上 右 侧 〉 标 注 在 分 支 上 ， 支 持 率 和 后 验 率 用 逗号 隔 开 。 滇 白 珠 复合 群 相关 变种 及 外 类 
群 系统 位 置 通过 不 同 颜 色 标识 : FX 为 外 类 群芳 香 白 珠 〈 蓝 色 )，WD 和 YHA EAA CLE), HXQ 和 
F HARAR CRE), WH 为 白果 白 珠 原 变 种 (黄色 )， 其 余 为 汗 白 珠 ( 黑 色 )。 分 支 A 包括 HXQ IF 居 
群 所 有 个 体 ;分 支 B1 包括 WH 居 群 所 有 个 体 , 4) 3c B2 包括 DYS M DYD 居 群 部 分 个 体 ,分支 B3 包括 FDQS 
和 FMHS 居 群 部 分 个 体 ， Arc B4 包括 云南 和 四 川 的 AN、DL、QJ、WD、EMS、XC、SB 居 群 的 全 部 或 部 
分 个 体 ; Ax Cl 包括 DWS 的 部 分 个 体 ， 分 支 C2 包括 FGT 和 FHJS 居 群 的 部 分 个 体 ， 分 支 C4 包括 SWT 
和 ZY 的 部 分 个 体 ， 剩 下 的 居 群 构成 一 个 大 的 分 支 C3。 系 统 树 下 部 展示 了 白果 白 珠 原 变 种 和 演 白 珠 复合 群 
3 个 国产 变种 的 植株 照片 。 


Fig. 1 Phylogenetic tree of the Gaultheria leucocarpa var. yunnanensis complex based on the Bayesian inference 


ug 


Ic 


(BD), with the population of G. fragrantissima as outgroup. The Bootstraps from the ML analysis (tree not shown) 
and the posterior probabilities are both shown on the branches. Numerals on the left side of branches indicate the 
bootstrap value (> 50% ) and those on the right side indicate the posterior probability value (> 5096), and the BP 
and PP values are separated with comma. The placements of populations of varieties within the complex and 
outgroup are indicated by branches in various color: FX (G. fragrantissima) in blue, WD and YHA (G. leucocarpa 
var. crenulata ) in red, HXQ and F (G. leucocarpa var. psilocarpa) in green, WH (G. leucocarpa var. leucocarpa) 


in yellow, and all the rest (G. leucocarpa var. yunnanensis) in black. Clade A includes all individuals of the HXQ 
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and F populations; clade B1 includes all individuals of the WH population; clade B2 includes partial individuals of 
the DYS and DYD populations; clade B3 includes partial individuals of the FDQS and FMHS populations; clade 
B4 includes all or partial individuals of the AN, DL, QJ, WD, EMS, XC, and SB populations; clade C1 includes 
partial individuals of the DWS population, clade C2 includes partial individuals of the FGT and FHJS populations, 
and clade C4 includes partial individuals of the SWT and ZY populations, and the remaining individuals form a 
large clade C3. Images beneath the phylogenetic tree show the close-up of plants of four varieties (three of them 


from China) within the G. /eucocarpa var. yunnanensis complex. 


Gaultheria fragrantissima 
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the clade from the Southeast China 
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图 2 汗 白 珠 复合 群 的 系统 发 育 网 络 图 ， 外 类 和 群 为 芳香 白 珠 。 大 于 或 等 于 50% 的 支持 率 显 示 在 分 支 上 。 主 要 
分 为 三 个 支 系 ， 即 : 横断 山脉 + 台湾 + 菲律宾 + 马来西亚 支 系 〈 红 色 圈 注 ， 右 侧 )、 华 东南 支 系 〈 蓝 色 圈 注 ， 
左 侧 )、 大 围 山 支 系 (紫色 圈 注 ， 中 央 )。 其 中 横断 山脉 分 支 居 群 包括 : 云南 的 AN、WD ERAR). 、QJ、 
SB、DL， 四 川 的 XC、EMS 居 群 ， 华 东南 支 系 包括 : 云南 的 PB、XS、ZY、TC、BS、MG、JD、YJ、WX、 
FY， 贵 州 ( 包 括 毛 溃 白 球 YHA 居 群 ，、 广 西 、 广 东 、 湖 南 、 湖 北 、 重 庆 、 江 西 、 福 建 的 所 有 居 群 ， 大 围 
山 支 系 仅 包括 云南 的 DWS 居 群 。 


Fig. 2 The Neighbor-Net Network of populations of the Gaultheria leucocarpa var. yunnanensis complex and 
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related taxa, with the population of G. fragrantissima as outgroup. Bootstrap values (2 50%) based on 1000 
replicates are indicated above branches. Three major clades are divided, i.e., the clade from Taiwan, the 
Philippines, Malaysia plus the Hengduan Mountains (colored in red, on the right), the clade from the 
Southeastern China (colored in blue, on the left), and the clade from Dawei Mountain with only one population 


DWS included (colored in purple, in the middle). Populations from the Hengduan Mountains clade include AN, 
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WD (G. leucocarpa var. crenulata), QJ, SB, and DL from Yunnan and XC and EMS from Sichuan. The clade from 
the Southeastern China is composed of populations from Yunnan (PB, XS, ZY, TC, BS, MG, JD, YJ, WX, and FY) 
and all the populations from Guizhou (include YHA of G. leucocarpa var. crenulata), Guangxi, Guangdong, 


Hunan, Hubei, Chongqing, Jiangxi, and Fujian provinces. 


3 讨论 


3.1 演 白 珠 复合 群 的 遗传 变异 

物种 的 形成 和 分 化 通常 源 于 群体 间 建 立 的 基因 流 屏障 ,主要 有 : 遗传 分 化 造成 的 生殖 隔 
离 、 地 质 及 气候 历史 变迁 造成 的 地 理 隔离 、 繁 育 系统 的 复杂 性 等 (Levin, 2003; Mallet, 2007; 
Rieseberg & Willis，2007)。 溃 白 珠 复合 群 分 布 广泛 ， 遗 传 差异 相对 形态 来 说 ， 与 地 理 分 布 
相关 性 更 大 。 从 地 形 上 来 看 ， 它 横 跨 了 二 、 三 级 阶梯 ， 从 横断 山脉 南端 一 直 治 长 江 流域 分 布 
到 东南 沿海 地 区 ， 海 拔 上 达 2 500 米 ， 下 至 390 米 。 野 外 调查 发 现 ， 该 植物 大 多 生长 在 空旷 
的 土 坡 或 岩 壁 上 , 呈 带 状 局 限 分 布 于 较 高 海拔 地 带 。 相 关 研 究 也 发 现 溃 白 珠 主 要 以 根茎 进行 
无 性 繁殖 ， 由 根 延伸 出 芽 形成 新 植株 〈 马 小 军 等 ，2001)， 且 马来西亚 的 白果 白 珠 原 变 种 为 
雄 营 先 熟 和 自 花 授 粉 (Sleumer，1967)。 所 以 滇 和 白 珠 复合 群星 现 出 的 遗传 分 化 可 能 是 地 形 的 
阻隔 导致 了 地 理 隔离 ， 处 于 横断 山脉 的 居 群 受 高 大 山脉 如 衣 牢 山 、 点 苍山 、 无 量 山 、 怒 山 、 
大 凉山 等 阻隔 ,而 处 于 东南 地 区 的 滇 白 珠 居 群 分 布 在 海拔 相对 较 低 的 大 起 伏 中 山 如 南 岭 、 大 
瑶 山 、 十 万 大 山 、 云 开 大 山 等 ， 地 势 较 横断 山脉 相对 为 平坦 和 连绵 〈 尤 联 元 等 ，2013 )， 沉 
白 珠 复 合群 可 能 因此 没有 产生 较 大 的 遗传 变异 。 


32 演 白 珠 复 合群 的 系统 演化 关系 
本 研究 采用 的 两 个 叶绿体 片段 rp133-psaJ 和 trnL-rpI32 在 滇 白 珠 复合 群 中 变异 位 点 较 
少 , 虽然 有 多 个 地 区 的 序列 没有 形成 分 支 , 但 解决 了 部 分 群体 的 系统 演化 关系 。 台湾 (HXQ) 
和 菲律宾 F) 的 秀 果 白 珠 居 群 得 到 了 分 子 系统 学 较 好 的 支持 (图 1 )， 从 遗传 上 支持 了 该 变 
种 的 成 立 ， 也 支持 了 子 房 壁 及 茹 果 不 被 毛 是 划分 该 变种 重要 的 共 衍 征 ( 分 类 性 状 )， 果 实 形 
e 态 特征 被 认为 在 白 珠 树 属 中 具有 重要 的 分 类 学 意义 (u etal, 20100 。 系 统 发 育 分 析 显 示 秃 
一 果 白 珠 是 汗 和 白 珠 复合 群 中 较 早 分 化 出 来 的 变种 , 但 和 具有 地 理 隔 离 的 横断 山脉 居 群 更 近 缘 的 
Oo 关系 并 没有 得 到 较 好 的 支持 (支持 率 低 )， 因 此 ， 与 横断 山脉 区 域 居 群 之 间 的 系统 发 育 关系 
还 需要 进一步 的 研究 。 在 1991 年 中 文 版 《中 国 植物 志 》 中 ， 屏 边 白 珠 因 其 叶片 形态 和 质地 
的 差异 被 作为 白果 白果 的 一 个 变种 处 理 ， 后 中 国 植物 志 英 文 版 (Fang & Stevens，2005) 将 
其 作为 滇 白 珠 的 异 名 处 理 ， 但 分 子 系统 发 育 显 示 屏 边 大 围 山 的 两 个 个 体 DWS-1 和 DWS-2 
与 其 他 滇 白 珠 居 群 分 属于 两 个 分 支 ， 并 且 位 于 较 基 部 而 较 早 地 分 化 出 来 〈 图 1); 从 网 状 进 
化 关系 看 ， 屏 边 大 围 山 的 居 群 被 单独 划分 为 一 文 〈 图 2)， 一 定 程度 上 支持 了 该 变种 遗传 学 
上 的 独立 性 。 
基于 传统 形态 分 类 ， 毛 滇 白 珠 因 其 枝 、 花 梗 密 被 均匀 分 布 的 硬 毛 ， 叶 边缘 被 缘 毛 而 被 单 
独 作为 一 个 变种 ， 区 别 于 汗 白 珠 (Fang & Stevens，2005)， 但 分 子 系统 发 育 显示 采 自 云南 武 
定 的 毛 滇 白 珠 和 采 自 云南 安宁 的 汗 白 珠 聚 为 一 支 , 同时 采 自 贵州 印 江 的 毛 汗 白 珠 和 同 地 区 的 
滇 白 珠 聚 为 一 支 〈 图 1)。 环 境 因素 对 植物 的 生长 发 育 会 产生 一 定 的 影响 ， 表 现 出 丰富 的 形 
态 、 和 生理 和 发 育 等 表 型 (Sultan，2000)。 根 据 我 们 野外 观察 发 现 生 活 在 背 阴 环境 的 汗 白 珠 
居 群 被 毛 情况 多 于 向 阳 环 境 ， 且 叶片 较 小 , 植株 较 纤细 。 本 研究 分 子 证 据 表 明 同 域 分 布 的 有 
毛 被 及 无 毛 被 汗 白 珠 在 遗传 上 没有 明显 差异 ， 说 明 被 毛 这 一 性 状 与 遗传 分 化 的 相关 性 不 大 。 
形态 与 遗传 的 不 一 致 性 在 许多 植物 中 也 有 报道 ， 例 如 : fe SEMESTRE, SEEMS 
Dasiphora floribunda (Pursh) Raf. 和 银 露 梅 D. glabra (G. Lodd.) Soj k 因为 花瓣 颜色 的 不 同 而 


被 分 为 两 个 物种 , 但 研究 发 现 基于 花 办 颜色 对 金 露 梅 属 植物 进行 物种 划分 界定 没有 得 到 遗传 
数据 的 文 持 CMa etal., 20140. Jib. AA ÆW E NE RIT EZB Haberlea 
ferdinandi-coburgii Urum 和 H. rhodopensis Friv 也 出 现 了 形态 变异 和 遗传 分 化 不 一 致 的 现象 

(Petrova et al，2014)。 对 于 汗 白 珠 复 合群 这 样 的 遗传 和 形态 不 一 致 性 ， 我 们 推测 可 能 的 原 
Hæ: C) 滇 白 珠 毛 被 的 生长 发 育 容 易 受到 环境 的 影响 , 这 种 变异 不 稳定 , 与 谱系 遗传 无 关 ， 
毛 滇 和 白 珠 只 能 作为 汗 白 珠 的 一 种 生态 类 型 ; (2) 叶绿体 片段 没有 足够 的 分 辨 率 且 仅 为 母系 遗 
传 ,无 法 检测 是 否 存在 杂交 和 网 状 进 化 ,需要 和 旬 选 合适 的 核 基 因 片 段 或 开展 简化 基因 组 分 析 
进行 更 深入 的 研究 。 


3.3 演 白 珠 复合 群 的 种 质 资源 保护 和 利用 

尽管 演 白 珠 复合 群 分 布 广泛 而 资源 丰富 ， 在 一 定时 期 内 可 以 承受 一 定 规模 的 开发 利用 。 
但 是 ， 该 植物 并 非 栖 生地 的 优势 种 ， 竞 争 性 差 、 生 长 相对 缓慢 (灌木 )、 生 物 量 有 限 、 气 候 
变化 敏感 性 高 ， 连 年 的 工业 化 采 控 会 导致 区 域 性 资源 很 快 枯竭 。 根据 野外 调查 发 现 ， 云 南 
部 、 东 北部 和 四 川 大 凉山 的 居 群 正 处 于 严重 的 人 为 干扰 中 , 该 区 域 的 演 白 珠 复 合群 经 历 了 当 
地 一 些 药材 公司 连续 性 整 株 ( 包 括 地 下 根部 ) 收购 ， 并 伴 有 民间 药 用 采 控 ， 这 些 居 群 的 生物 
量 正 面临 不 同 程度 的 减少 , 对 该 复合 群 进行 保护 生物 学 的 研究 极其 重要 , 需要 在 本 研究 的 基 
础 上 , 进行 居 群 式 取样 开展 生态 地 理学 及 谱系 地 理学 研究 , 更 深层 次 阐明 遗传 变异 格局 的 形 
成 机 制 ， 确 定 保护 单元 ， 才 能 制订 有 效 的 保护 策略 。 
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